Subcritical finite-amplitude solutions in plane Couette flow of visco-elastic fluids 
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Plane Couette flow of visco-elastic fluids is shown to exhibit a purely elastic subcritical instability 
in spite of being linearly stable. The mechanism of this instability is proposed and the nonlinear 
stability analysis of plane Couette flow of the Upper-Convected Maxwell fluid is presented. It is 
found that above the critical Weissenberg number, a small finite-size perturbation is sufficient to 
create a secondary flow, and the threshold value for the amplitude of the perturbation decreases as 
the Weissenberg number increases. The results suggest a scenario for weakly turbulent visco-elastic 
flow which is similar to the one for Newtonian fluids as a function of Reynolds number. 
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Recently, it has been discovered that flows of visco- 
elastic fluids like polymer solutions and melts can lose 
their stability and become turbulent at very low Reynolds 
numbers 0,12, In contrast to the Newtonian turbu- 

lence, where inertia plays a destabilizing role, this elastic 
turbulence or turbulence without inertia Q arises when 
the flow-induced stretching of the polymers makes the 
elastic stresses in the fluid become large and anisotropic. 
It is a challenge to explain how the transition to this new 
type of turbulence can occur. In this Letter we show that 
visco-elastic plane Couette flow (PCF) exhibits a subcrit- 
ical instability to finite amplitude waves and argue for a 
scenario for weak elastic turbulence analogous to that for 
weak turbulence in PCF of Newtonian fluids. 

The mechanism of linear elastic instability was identi- 
fied for flows with curved stream-lines 0, |(| . One of the 
classical examples of such a flow is realized in Taylor- 
Couette cell where fluid fills the gap between two coax- 
ial cylinders made to rotate with respect to each other. 
In the laminar state, the fluid moves around the cylin- 
der axis and the elastic or hoop stresses act on poly- 
mer molecules stretching them along the circular stream- 
lines and exerting extra pressure towards the inner cylin- 
der. When these stresses overcome viscous friction, the 
laminar state becomes linearly unstable - any infinitesi- 
mal perturbation will push a polymer from the circular 
stream-lines and create a secondary flow, the so-called 
Taylor vortices. Pakdel and McKinley generalized this 
mechanism to arbitrary flows Q and proposed that there 
exists a universal relation between the properties of the 
fluid and the flow geometry which determines the condi- 
tions of the linear instability. One of the dimensionless 
parameters in their argument is the so-called Weissenberg 
number Wi = A 7, where A is the elastic relaxation time 
of the fluid, and the shear rate 7 = dv/dx gives the 
relative velocity of two fluid layers moving with respect 
to each other. Pakdel and McKinley argued that 
the critical Weissenberg number is related to the char- 
acteristic curvature of the flow stream-lines and that the 
linear instability disappears when the curvature goes to 
zero. The known results on the visco-elastic instabilities 



in Taylor-Couette yj, cone-and-plate and parallel plate 
, Dean and Taylor-Dean [|| flows are in agreement with 
this curved stream-lines - linear instability paradigm. 

The linear stability of parallel visco-elastic shear flows 
has been investigated in detail. For essentially all stud- 
ied visco-elastic models, laminar PCF is linearly stable 
HEI (note the exception 0). In the case of the pipe 
flow, the linear stability was demonstrated numerically 
by Ho and Denn for any value of the Weissenberg 
and Reynolds numbers. Therefore, it has become com- 
mon knowledge that the parallel shear flows of fluids 
obeying simple visco-elastic models (UCM, Oldroyd-B, 
etc.) are linearly stable, in agreement with the curved 
stream-lines - linear instability paradigm. Clearly, if an 
instability does occur in practice, it has to be nonlinear. 

There is no general agreement, however, on whether 
parallel shear flows like PCF or pipe do in fact become 
unstable. At the moment, there has been no experiment 
that would clearly establish the presence or absence of 
a bulk hydrodynamic instability in parallel visco-elastic 
shear flows. One of a few indirect indications that a bulk 
instability might occur in the pipe flow comes from the 
famous melt-fracture problem |l2| . which arises in ex- 
trusion of a dense polymer solution or melt through a 
thin capillary. There, when the extrusion rate exceeds 
some critical value, the surface of the extrudate becomes 
distorted and the extrudate might even break, giving 
the name to the phenomenon. It is possible that this 
is a manifestation of an instability taking place inside 
the capillary, though other mechanisms (such as stick- 
slip, influence of the inlet, etc.) have been proposed [l^ . 
Recently, we presented arguments for the bulk instabil- 
ity being related to the melt-fracture phenomenon |l3| . 
but the issue stays highly controversial. There is also 
some evidence for nonlinear parallel shear flow instabil- 
ities from numerical simulations of visco-elastic hydro- 
dynamic equations 14] . Partly because the numerical 
schemes used to solve these equations are known to break 
down when elastic stresses become large (Wi > 1) — the 
so-called large Weissenberg number problem — it is 
open to debate whether an observed phenomenon is due 
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to a numerical or a true physical instability. 

In this Letter we show explicitly that a nonlinear in- 
stability mechanism does exist in agreement with the fol- 
lowing argument. The laminar velocity profiles of the 
parallel shear flows have straight stream-lines, and, there- 
fore, their linear stability is in agreement with the curved 
stream-lines - linear instability paradigm. The linear the- 
ory predicts that a small perturbation superimposed on 
top of the laminar flow will decay in time with the de- 
cay rate depending on the Weissenberg number. When 
Wi becomes larger than one, the decay time becomes 
comparable with the elastic relaxation time A, and the 
perturbation becomes long-living. Thus, on short time- 
scales, the superposition of the laminar flow and the 
slowly-decaying perturbation can be viewed as a new ba- 
sis profile with curved stream-lines. Applying the same 
curved stream-lines - linear instability paradigm to the 
perturbed streamlines, we conclude that this new flow 
can become linearly unstable. The instability occurs 
via a perturbation on top of a perturbation, and thus 
is nonlinear. Since the initial perturbation has to be 
strong enough to become unstable, there exists a finite- 
amplitude threshold for the transition, which becomes 
smaller as the Weissenberg number increases. This sce- 
nario resembles transition to turbulence in parallel shear 
flows of Newtonian fluids. There as well, one encoun- 
ters the absence of the linear instability, and a subcrit- 
ical transition with the threshold going down with the 
Reynolds number 0, 0] . 

Our explicit results are for the nonlinear stability anal- 
ysis of PCF. We consider the so-called UCM fluid 
confined in the y-direction in-between two plates (y = 
id), which move with constant velocity vq in the oppo- 
site directions along the cc-axis. The hydrodynamic equa- 
tions, consisting of the equations for momentum balance, 
the UCM model, and incompressibility, read (lif 
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= -[(Vv) + (Vv)t], (2) 
V-v = 0, (3) 

where the Weissenberg number Wi=Xvo/d, the Reynolds 
number Re = pv^/(rjd), p is the density and r\ is the 
viscosity of the fluid, r and v are the stress tensor and the 
velocity, respectively; d is used as the unit of length, d/vo 
as the unit of time, and the stress tensor is scaled with 
ijvo/d; (...y denotes the transposed matrix. As usual, 
we split these variables in two parts: the laminar values 
{T~iam, y^x} and the deviation describing the disturbance 
{t',v'}. It is useful to organize all hydrodynamic fields 
of the disturbance in one vector V = {t[ , v[,p}^ . Then, 



the equations illl.'SI) can formally be re-written as 
d 

C V + — {Re v', Wi t', 0} t = N {V, V) , 
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where the l.h.s. represents the linear terms in eQS. (|llH| l. 
and the r.h.s. the quadratic nonlinearity. 

The first step of our analysis is to determine the eigen- 
values A and the eigenfunctions Vq of the linear operator 
L. Gorodtsov and Leonov and Reynardy et al. 
have shown that the eigenfunctions in the form 



V (x,y,z) = V Q (y)e^ kx+ ' 1 ^ + c.c. 



(5) 



have two types of physical eigenvalues: a pair of 
complex-conjugated "elastic" eigenvalues (Gorodtsov- 
Leonov modes), and an infinite discrete set of "inertial" 
ones. (There also exists a continuous spectrum of eigen- 
values which were shown to be unphysical [l9| and which, 
therefore, will be discarded.) Let us for the moment fo- 
cus on the elastic mode, and suppose we choose the initial 
disturbance to be in the form of the elastic eigenmode 

V{x,y,z,t) = <P(t)V (GL) (y)e^ kx+qz '> + c.c. (6) 

where $(t) is a complex amplitude (Our normalization 
is such that when the amplitude <I>(i) = 1, the strength 
of the shear rate created by the perturbation equals that 
of the laminar flow) . In order to investigate how the am- 
plitude $(t) will change in time depending on Wi, k and 
q, we derive the equation governing the time evolution of 
<&(t), or the amplitude equation. The standard technique 
used to derive the amplitude equations relies on the pres- 
ence of linear instability (occurring at, say, Win n ), and 
uses the distance to the instability (Wi — Wi\i n )/Wiu n as 
a small parameter. Then, the nonlinear evolution of <I>(i) 
near Wi^ can be deduced by the method of multiple 
scales 20] . The parallel shear flows, however, are lin- 
early stable (Winn — > oo, effectively), and we have to use 
other methods. Instead, assuming $(t) to be small, we 
substitute Eq. JHJ into Eq. collect the terms propor- 
tional to exp [i(kx + qz)], and using the corresponding 
eigenmode of the adjoint operator C' , we project these 
terms on the original form Eq. © . We then find the time 
derivative of $(£) as a series in $(i) 
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where A' GL - ) is the eigenvalue of the elastic (Gorodtsov- 
Lconov) mode, and the nonlinear coefficients C"s are 
explicit functions of Wi, k and q. Eq. @ has solu- 
tion $(t) = |$|e lf2 ' which results in traveling waves in 
Eq. jnj). For small amplitudes, it reproduces the linear 
decay $ ~ exp [\( GL h] , Re (\( GL 1) < 0. The instability 
threshold is determined by finding the steady-state solu- 
tions of eq.0. The main problem in dealing with the 



FIG. 1: Steady-state amplitude <E> for fc = 1 and q = 1: 
dashed, dot-dashed, dotted and solid lines show the solution 
to Re = up to the 5 th , 7 th , 9 th and ll" 1 order. The 
inset illustrates qualitatively that the series converges only 
for small enough |$|. The ratio Re/Wi = 10 -3 was kept 
constant. 
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FIG. 2: The saddle-node Weissenber number as a function 
of the wave-vectors. The dots are the calculated values, the 
lines serve as a guide to the eye. 



series like eq. (JJJ is that it is not known a priori whether 
it converges. In order to check the convergence of the 
series upon inclusion of higher-order terms, we solve the 
equation Re (=£) — in successive orders. 

In FigQ] we show the solution to these equations for 
k = 1 and q = 1. The most important feature of these 
curves is that they show existence of a subcritical insta- 
bility for Weissenberg numbers larger than the saddle- 
node value Wi sn indicated in the figure and obtained 
by extrapolating the curves. As the arrows indicate, for 
Wi > Wi sn the lower branch of the curves denotes the 
critical amplitude — amplitudes larger than this value 
will grow in time. Note that the instability threshold 
is small (consistent with our assumption |$| < 1), and 
goes down as Wi increases. The inclusion of higher-order 
terms causes the whole curve to shift to the right, though 
the shift becomes roughly two times smaller with every 
coefficient included, indicating convergence. 

While the lower branch of each curves gives the min- 
imal amplitude of the disturbance sufficient to destabi- 
lize the laminar flow, the upper branch determines the 
saturated value of $ after the transition. Surprisingly, 
it diverges in the vicinity of the saddle-node where the 
highest coefficient in the expansion changes sign. There 
could be several reasons for that. First, it may indicate 
that the nonlinear state in the form of Eq. © is unstable 
and will undergo a transition to another coherent state 
or to turbulence. Second, although we consider this un- 
likely, it may be that the UCM model cannot capture 
this state. Finally, it may be a convergence problem. As 
indicated in the inset of Fig^ the upper branch may lie 
beyond the radius of convergence of Q . When we apply 
the method to the subcritical Swift- Hohenberg equation 
|2lj . we also find such behavior. 

In FigEl we plot the lowest Weissenberg number for 
which the nonlinear instability is possible, or the position 



of the saddle- node Wi STl , as a function of the wave- vectors 
k and q. It clearly shows that the saddle-node position 
is only a weak function of the wave-vectors, and a large 
number of modes with different fc's and q's is nonlincarly 
unstable for given Wi > 2.1. Even if each individual 
mode saturates at a given value of <&, the superposition of 
a large number of such modes is likely to become chaotic, 
and we expect the elastic turbulence to set close to or 
even at the instability. 

The existence of the subcritical instability in PCF can 
also be deduced from the well-studied Taylor-Couette 
flow 0, [|| [2^. There, the subcritical character of the 
transition was both observed experimentally |2^ | and cal- 
culated analytically [24| . We have performed the nonlin- 
ear stability analysis for the purely elastic eigenvalue of 
the non-axisymmetric mode used by Sureshkumar et al. 
[2^ |. which we have traced to the case of the counter- 
rotating cylinders. We have found that this mode be- 
comes subcritically unstable far below the linear insta- 
bility threshold. Moreover, as the curvature of the cylin- 
ders goes to zero, the linear instability disappears at in- 
finity, while the threshold of the subcritical transition 
goes smoothly to the PCF one calculated above. This 
is somewhat similar to the Nagata's construction for the 
Newtonian fluids [2^. The full report on these results 
will be published elsewhere. 

Now we turn to the discussion of the " inertial" eigen- 
modes. Replacing A' Gi - ) and Vq GL ^ in Eq. © with one 
of the inertial eigenvalues and the corresponding eigen- 
mode, and repeating the same calculation for the coef- 
ficients in the amplitude equation @, we find that this 
type of disturbances is nonlinear ly stable for Wi ~ Wi sn . 
We have checked this result for the three eigenvalues with 
the smallest imaginary parts. The structure of the eigen- 
spectrum also helps to rule out the possibility that the 
disturbances in the form of the inertial eigenmodes can 
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become unstable via a nonlinear interaction with the elas- 
tic mode. To the lowest order, the interaction between 
the elastic amplitude <I> and an inertial amplitude 'J 7 is 
described by adding terms like Cg**^$|^| 2 to the ampli- 
tude equation (JJJ, and similar terms to the equation for 
Since the imaginary part of the elastic eigenvalue is 
0(1), while the imaginary part of the inertial eigenvalues 
is 0(1/V Re) ^> 1 for small Re, the coefficient C3 , 
describing the overlap between two modes, must be very 
small, so the interaction terms can be neglected. 

The above results give very strong indications for the 
existence of a branch of nonlinear finite-amplitude so- 
lutions which renders visco-elastic PCF flow nonlinearly 
unstable for Wi > 2.1. We believe these solutions do play 
an important role in organizing the dynamics of visco- 
elastic PCF. First of all, it is intriguing to note that in the 
direct numerical simulations of PCF of Atalik and Keun- 
ings [l4j , numerical instabilities were found for Wi > 2 — 
could this be due to the occurrence of these modes? Sec- 
ondly, we have found that the branch of solutions which 
we have established here have an analog in the pipe flow 
where the transition is found to be at Wi sn w 5 (results 
will be published elsewhere, see also 13] ). Thirdly, and 
most importantly, the similarity of our results with what 
has been found for PCF of Newtonian fluids lead us to 
suspect that these solutions could be the building blocks 
of weak visco-elastic turbulence: The visco-elastic peri- 
odic finite-amplitude solutions found here resemble the 
periodic finite-amplitude solutions in Newtonian fluids. 
There, they are known to be exact but unstable solu- 
tions of the Navier-Stokes equation [2^ and take part in 
the self-sustaining cycle — a periodic sequence of insta- 
bilities in which streaks, stream-wise vortices and rolls 
are continuously destroyed and regenerated [2H. EtI ] . It 
is tantalizing to speculate that a similar cycle can be 
proposed to sustain the weak elastic turbulence in PCF 
and Poiseuille flow. The first step in this direction was 
made in |28j where the influence of a minute amount of 
polymer on the weak Newtonian turbulence was studied. 
Together with our findings on the finite-amplitude solu- 
tions this gives hope that a visco-elastic version of the 
self-sustaining cycle can indeed be formulated. 

In conclusion, we have presented evidence for nonlin- 
ear instability in visco-elastic PCF. Together with the 
vast experimental and numerical evidence in various flow 
geometries 0, S E this makes us believe 
that the nonlinear subcritical instability is an inherent 
feature of the visco-elastic parallel shear flows and that 
the finite-amplitude solutions may organize dynamics of 
weak elastic turbulence. 



[1] R. G. Larson, E. S. G. Shaqfeh, and S. J. Muller, J. Fluid 
Mech. 218, 573 (1990). 



[2] G. H. McKinley, J. A. Byars, R. A. Brown, and R. C. 
Armstrong, J. Non-Newtonian Fluid Mech. 40, 201 

(1991) . 

[3] A. Groisman and V. Steinberg, Nature 405, 53 (2000). 
[4] A. Groisman and V. Steinberg, New J. Phys. 6, 29 (2004). 
[5] R. Larson, Nature 405, 27 (2000). 

[6] Y. L. Joo and E. S. G. Shaqfeh, Phys. Fluids A 4, 524 

(1992) . 

[7] P. Pakdel and G. H. McKinley, Phys. Rev. Lett. 77, 2459 

(1996) . 

[8] V. A. Gorodtsov and A. I. Leonov, J. Appl. Math. Mech. 
31, 310 (1967). 

[9] M. Renardy and Y. Renardy, J. Non-Newtonian Fluid 
Mech. 22, 23 (1986); M. Renardy, Eur. J. Mech. B 11, 
511 (1992); H. J. Wilson, M. Renardy, and Y. Renardy, 
J. Non-Newtonian Fluid Mech. 80, 251 (1999). 
[10] A. M. Grillet, A. C. B. Bogaerds, G. W. M. Peters, and 
F. P. T. Baaijens, J. Non-Newtonian Fluid Mech. 103, 
221 (2002). 

[11] T. C. Ho and M. M. Denn, J. Non-Newtonian Fluid 
Mech. 3, 179 (1977/1978). 

[12] C. J. S. Petrie and M. M. Denn, AICHE J. 22, 209 (1976); 
M. M. Denn, Annu. Rev. Fluid Mech. 22, 13 (1990); 
M. M. Denn, Annu. Rev. Fluid Mech. 33, 265 (2001). 

[13] B. Meulenbroek, C. Storm, V. Bertola, C. Wagner, 
D. Bonn, and W. van Saarloos, Phys. Rev. Lett. 90, 
024502 (2003); V. Bertola, B. Meulenbroek, C. Wag- 
ner, C. Storm, A. N. Morozov, W. van Saarloos, and 
D. Bonn, Phys. Rev. Lett. 90, 114502 (2003); B. Meu- 
lenbroek, C. Storm, A. N. Morozov, and W. van Saarloos, 
J. Non-Newtonian Fluid Mech. 116, 235 (2004). 

[14] K. Atalik and R. Keunings, J. Non-Newtonian Fluid 
Mech. 102, 299 (2002). 

[15] R. G. Owens and T. N. Phillips, Computational Rheology 
(Imperial College Press, 2002). 

[16] P. J. Schmid and D. S. Henningson, Stability and transi- 
tion in shear flows ( Springer- Verlag, New York, 2001). 

[17] B. Hof, A. Juel, and T. Mullin, Phys. Rev. Lett. 91, 
244502 (2003). 

[18] R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics 
of polymeric liquids, vol. 1 (John Wiley & Sons, Inc., 
1987), 2nd ed. 

[19] M. D. Graham, J. Fluid Mech. 360, 341 (1998). 

[20] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

[21] A. N. Morozov, P. Becherer, and W. van Saarloos, in 
preparation. 

[22] M. Avgousti and A. N. Beris, J. Non-Newtonian Fluid 

Mech. 50, 225 (1993). 
[23] A. Groisman and V. Steinberg, Phys. Rev. Lett. 78, 1460 

(1997) . 

[24] R. Sureshkumar, A. N. Beris, and M. Avgousti, Proc. R. 
Soc. Lond. A 447, 135 (1994). 

[25] M. Nagata, J. Fluid Mech. 217, 519 (1990). 

[26] F. Waleffe, Phys. Fluids 9, 883 (1997). 

[27] B. Hof, C. van Doorne, J. Westerweel, F. Nieuwstadt, 
H. Faisst, B. Eckhardt, H. Wedin, R. Kerswell, and 
F. Waleffe, Science 305, 851 (2004). 

[28] P. A. Stone, F. Waleffe, and M. D. Graham, Phys. Rev. 
Lett. 89, 208301 (2002); P. A. Stone, A. Roy, R. G. Lar- 
son, F. Waleffe, and M. D. Graham, Phys. Fluids 16, 
3470 (2004). 

[29] K. A. Kumar and M. Graham, Phys. Rev. Lett. 85, 4056 
(2000). 



